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Abstract. Observations carried out in the frame of the ISTC project 3547 [1] 
were used to parameterize and validate the model of radionuclides transport 
from Siberian Chemical Combine (SCC) by the Tom and Ob rivers and also to 
assess discharges of radioactive substances from SCC to the Tom River and to 
estimate possible contamination of the rivers in case of accidents. 


1. INTRODUCTION 
The SCC was founded about 60 years ago and is the biggest Nuclear Fuel Cycle 
plant in Russia. The main sources of radioactive contamination of the Tom River 
were its industrial reactors. The last of the reactors was closed in 2008. 
Contaminated water discharges from SCC were carried out through the 
Romashka River to the Chernilschikovskaya channel - the right distributary of 
the Tom River [2-8]. Observations of specific activities of radioactive substances 
in the Tom and Ob rivers were used to assess discharges of radioactive 
substances from the SCC. The assessment was carried out with the use of the 
model based on the two-dimensional equation of advection and dispersion. It 
takes into account exchange of radionuclides between water column (solute, 
suspended particles) and bottom sediments [9]. 
In current work the following additional assumptions were made: 
e Influence of longitudinal dispersion is negligible in comparison with 
advection. Only lateral dispersion is taken into account; 
e Characteristics of a river channel are constant (see. Table 3). 


The system of differential equations describing migration of 
radionuclides is as follows: 
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* The current research is a part of the ISTC project 3547 


Where: ‘t’ is time, s; ‘F’ describes sources of a radionuclide, Bq/(m**s); ‘y’ 
is a distance from the right bank, m; ‘x’ is a distance downstream from the 
source of contamination, m; ‘Cw’ and ‘C,’ are the activities of a radionuclide per 
unit volume of water and bottom sediments respectively, Bq/m’; ‘A’ is the 
radioactive decay constant, st; ‘H’ is the average depth of a water object, m; ‘u’ 
is the flow velocity, m/s; ‘h’ is the thickness of the upper (effective) layer of 
bottom sediments, m; ‘B’ is the coefficient of diffusive mass exchange of a 
radionuclide between water column and upper layer of bottom sediments, m/s; 
‘y’ is the coefficient of diffusive mass exchange of a radionuclide between upper 
and lower layers of bottom sediments, m/s; ‘U’ is the effective sedimentation 
rate of suspended particles, m/s; Y is the intensity of resuspension, m/s; ‘opw’ 
and ‘ap’ are fractions of a radionuclide dissolved in water column and in upper 
layer of bottom sediments respectively, dimensionless; ‘atw’ and ‘am’ are 
fractions of a radionuclide sorbed on particles in water column and in upper 


layer of bottom sediments respectively, dimensionless; E, - the coefficient of 
lateral turbulent dispersion, m?/s. 


Values of Qm, Aq, pw Ap, can be determined with the use of well- 
known dependencies: 
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where ‘&, is the partitioning coefficient of a radionuclide between water and 


suspeneded matter, m°/kg; ‘ky - the partitioning coefficient of a radionuclide 


between poral water and solid phase, m°%/kg; ‘Sı’ - the concentration of 
suspended matter in water, kg/m; ‘m’ - the air-dry weight of unit volume of 


bottom sediments, gin. 
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Figure 1. The scheme of sampling regions (see Table 1). 


With the use of constants ‘Ay’, ‘Ao’, ‘Mz and ‘Az’ (1) can be written as 
follows: 
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Hydrological data about the Chernilschikovskaya channel are scarce. 
That’s why the width of the channel was estimated on the base of satellite 
photos (www.google.com) while flow velocity on the base of expression: 


u=C,\ RI, 
(5) 

Where ‘R’ is the hydraulic radius, which is approximately equal to 
average depth, m; ‘J’ is slope of a river; ‘Csa’ is the Shesi’s coefficient, which 
describes resistance of the friction between water flow and bottom sediments. 
The coefficient can be estimated with the use of the Manning’s expression. 


The empirical Karaushev’s expression was used to estimate ‘ poe 
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Where ‘B’ is the width of a river, m. 
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2. ASSESSMENT OF RADIOACTIVE DISCHARGES TO THE TOM RIVER 
DURING OPERATION OF THE INDUSTRIAL REACTORS 


Let us assume that discharges were steady and lasted for a long time. Let us 
also assume that the channel data and hydrological characteristics of the rivers 


OG... oC, 
were constant. Then time derivatives in (1) can be set to zero ae z =0). 
O 


Thus, for a steady-state problem, solution of (1) is as follows: 
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Here’ Cadorna - background activity of a radionuclide per unit 


(8) 


volume of water, Bq/m®; ‘A’ - intensity of inflow of a radionuclide to the Tom 
liña ) 

As ` 
Assessment of radioactive substances inflow intensity to the Tom River 


was carried out on the base of observed activities of the radionuclides in the 
rivers, with the use of the expression (9), derived from (7): 
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Here x, yi are coordinates of an observation spot, m; ‘C „ is the 


River, Bq/s; A=(1,- 


(9) 


observed activity, Bq/m?. The observed activities [1] one can find in Table 1. The 
regions of sampling are shown on Figure 1. 


Table 1. The observed activities that were used for the 
assessments. 


Sampling spot Date of | '8’Cs, Bq/m? Sr, 239,240 P: | 23ND, 
sampli | Suspend | filtrate | Bq/m? | u Bq/m’ 
ng ed filtrate | mBq/ | suspend 

particles mê ed 
filtrate | particle 
s 

Region 1 (r. 27.05. | 0.150. | 0.100. | 4.60. 

Tom, backgroun 08 02 02 7 

d) 24.06. | 0.080. 3.6+0. 

08 01 7 

25.07. 3,2+0, 

08 6 

01.09. | 0.06+0. 4,1+1, | 20+10 
08 02 0 

15.09. 3,8+1, 

08 1 

22.10 2,50, 

08 7 

Region 2 (r.|r. | 30.05 1.1340. | 0.240. | 8.541. | 60413 | 10645 

Tom, 8 km |b. | 08 04 03 0 

downstream 26.06. | 0.4740. | 0.540. | 11.1+ 

from the mouth 08 03 08 1.3 

of Romashka 01.08. | 0.340. | 0.150. | 12.1+ 

River) 08 03 08 1.5 
31.08 0.74+0 8,342, | 45+11 
08 04 1 

2.4+0.1 
29.09. 7,1+1, 
08 8 
m. | 31.08. | 0.30+0 


08 03 
Region 3 (r./r. | 29.05 1.14+0. | 0.15+0. | 6.241. | 55412 | 43+3 
Tom, 20 km |b. | 08 04 06 0 
downstream 27.06. 8.7+1 
from the mouth 08 1 
of Romashka 29.07 0.61+0. | 0.78+0. | 10.5 
River) 08 04 07 1.4 

29.08 0.70+0 0.670. | 6.2+1. | 41+9 

08 06 09 0 

24.09 5,7+1, 

08 5 

l.b | 29.08 0.12+0 
; 09 05 

Region 4 (r.J/r. |27- 0.360 2947 
Rom, 30 km | b. | 28.08 03 
downstream m. | 08 0.44+0 
from the mouth 05 
of Romashka | ].b 0.1640 
River) 02 
Region 5 (r. Ob, | r. | 24- 0.34+0 2845 
16 km b. | 25.08. | 04 
downstream m. | 08 0.18+0 
from the mouth 02 
of the Tom] 1.b 0.07+0 


Quantity of 72°Np in discharges of direct-flow reactors is always sufficient 
for reliable measurements. Data on actual discharges of this radioactive 
substance by SCC are available (see Table 2). That’s why observed data on 7°°Np 
were used for parameterization and validation of the model of the considered 
part of the Tom River. 
Table 2. - Intensity of discharges of 72°Np to the Tom River, Bq/year [2-8]. 


Actual discharge Permitted 
2002 2003 2004 2005 2006 2007 2008* discharge 
8.14-10 | 6.23-10 | 7.51-10 | 13.0-10 | 14.6-10 | 12.6-10 | 7.15-10 | 1.48-10"° 
12 12 12 12 12 12 12 


* - operation of direct-flow reactors were terminated 5.07.2008 


Assessment of the intensity of discharges carried out in accordance with 
(9) on the base of 7°°Np observed activity in water - (13.6+0.3):10” is in 
satisfactory agreement with data in Table 2. It enables one to use estimated 
values of characteristics of the Chernilschikovskaya channel for assessment of 
other radionuclides inflow intensity. In Tables 3 and 4 one can find values of 
parameters used for the assessments. In Table 5 one can find intensity of 
discharges assessed on the base of observed activities. 


Table 3. Channel data and hydrological characteristics used for 


modeling. 
Main channel of the Chernilschikovskaya Ob 
Tom River channel River 
B, m 700.0 450 1000 
H, m 4.0 4.0 6 


h,m 0.1 0.1 0.1 
Q, m°/s 980 630 3000 
Ey, m*/s 0.49 0.27 0.99 

Sı, 0.03 0.03 0.1 
kg/m? 

By, 1.9-10° 1.9:10° 1.9:10° 

8 
m/s 
U, m/s 1.0-10* 1.0-104 1.0:10° 
4 
Y, m/s 2.7-10° 2.7-108 9.1-10 
8 
m, 1100 1100 1100 
kg/m? 


Table 4. Estimated values of partitioning coefficients of the 
radionuclides and their background activity. 


Cs Sr 239.240 Pu 239Np 
k, 100 2 30 2.0 
k; 10 0.1 0.5 0.1 
Background activity in water, 0.1* 3.63** 0.02** 0 


Bq/m? 


* - on suspended particles, ** - in filtrate 


Table 5. Assessment of long-term discharges of the radionuclides to 
the Tom River, Bq/year. 


137Cs 


9Sr 239,240 Pu 


(2.2+1,.7)-10¥ (4.8+0.9)-10'° (6.5+1.9)-10° 


3. ESTIMATION OF POSSIBLE ACCIDENTAL CONTAMINATION OF THE 
TOM RIVER 
A scenario of hypothetical accident was studied in the frame of ISTC project 
3547. The scenario includes an accident followed by transport of radionuclides 
in the atmosphere and radioactive fall-out on the surface and the catchment 
area of the Tom River. A similar accident took place there on 06.04.1993 [10]. 
The most adverse fall-out conditions were assumed for assessment of the river 
contamination by long-living °°Sr and '%’Cs. As no detailed information about 
possible source of contamination was available, the modeling estimations were 
done for “unit” amount of radioactivity (1 TBq of each of the radionuclides). 

Modelling was carried out with the use of computer model Cassandra 
[11]. It enables computations in accordance with (1). A reach part of the Tom 
River with silty bed sediments and pronounced sedimentation, located 26-28 km 
downstream from the confluence with the Romashka River was taken into 
account. The reach reduces peak of specific activities of radionuclides 
(especially of '°’Cs) in the zone of contaminated water, moving downstream the 
river system. The reason is sedimentation of suspended particles, containing 
substantial fraction of '*’Cs, while the contaminated water is passing the reach. 
On the other hand afterward the silty bed sediments will become the source of 
secondary contamination of the water. 

On the Figure 2 one can see distribution of the radionuclides activity 
along the Tom River in 10 and 20 hours after the fall-out. The fall-out on the 
river surface was assumed to be instantaneous and to affect area of 2 km along 


the channel. Modelling results showed that after passing the reach mentioned 
above, peak activities of Cs reduced from 93 to 19 kBq/m? and °°Sr from 160 
to 108 kBq/m?. Maximum activity of '°’Cs (7 kBq/m®) and °°Sr (85 kBq/m?) in the 
mouth of the Tom River is to be observed in 32 hours after the fall-out. In 40 
days after the hypothetical accident, specific activity of '’Cs in water of the 
reach part of the Tom River is to be 20-40 Bq/m°. Specific activity of °°Sr is to be 


7-15 Bq/m? there. Maximum assessed specific activity of the radionuclides in the 
Ob River one can see in Table 6. 
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Figure 2. Distribution of radioactivity along the Tom River in 10 and 20 hours after the 
hypothetical fallout. 


Table 6. Maximum exceeding of background activity in the Ob River, Bq/m?. 


Right bank Left bank 
137Cs Sr 137Cs 9Sr 
16 km downstream from the mouth of the 5000 88000 0 0 
Tom River 
66 km downstream from the mouth of the 250 12600 66 1.5 
Tom River 
130 km downstream from the mouth of the 46 6200 370 3 
Tom River 


4. CONCLUSION 


The study enabled one to use observed data for assessments of long-living 
radionuclides discharges to the Tom River that took place during operation of 
the industrial reactors on the SCC. The hypothetical accident was studied. Its 
scenario included transport of Sr and '*’Cs in the atmosphere and their fall-out 
on the surface of the Tom River near the SCC. Modelling was carried out with 


the use of the computer model Cassandra. It provided estimation of spatial and 
temporal distribution of the radionuclides in the Tom and Ob rivers in case of 
fall-out of 1 TBq of each of the radionuclides. 
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